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D I Abstract: We checked Eling-Oz formula [1] for the bulk viscosity of the holographic fluid 

dual to the exotic black holes [2]. Initially, Eling and Oz argued that the formula is valid 
in the high temperature and adiabatic limits. In [3] the validity of the formula for J\f = 2* 
plasma and cascading gauge theory was pushed forward for arbitrary temperatures. We 



\f^ ' successfully verified the formula with the computations of the bulk viscosity [4, 5] for a 

OO ■ wide range of the temperatures. Moreover, it correctly reproduces the critical behavior in 

-vj . the vicinity of the critical point, where the bulk-to-shear viscosity diverges. 
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1 Introduction and Summary 

A new formula for bulk-to-shear viscosity of strongly coupled gauge theory plasma was 
proposed by Eling and Oz (EO) [1] . The wide class of the gauge theories are dual (in the 
context of AdS/CFT [6]) to the following {d+ 1) gravitational action 

S=-^lV^d''-^'x(R-^Y.(d<P,f-Vi<pA +S,auge. (LI) 

They used the null focusing (Raychaudhuri) equation describing the evolution of the horizon 
entropy, which is equivalent to the viscous fluid entropy balance law. In the absence of 
chemical potentials for the conserved charges, the formula for the bulk viscosity of the 
plasma dual to (1-1) takes the following form 
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where cj)^ are the scalar field values evaluated at the horizon of the black brain, T is the 
temperature of plasma dual to the black brane, Cs is the speed of sound waves in plasma. 
In the same paper the EO formula was verified for a large number of gauge theories dual 
to string theory at high temperature limit [7-10] and some phenomenological models of 
gauge/gravity correspondence [11, 12]. 

The expression for the bulk-to-shear viscosity employs the values of scalar fields only 
at the horizon. It is an intriguing result as the bulk viscosity in general depends on the 
energy scale; the boundary data is essential to capture microscopic scales of the theory. 
That is in contrast with the universality of the shear viscosity calculations [13-15]. In [3] 
the validity of (1.2) was extended for cascading gauge plasma [16, 17] and M = 2* gauge 
theory plasma [18-20] for all the temperatures. 

The correctness of the EO formula for the phenomenological models of gauge/gravity 
correspondence was also verified in [21]. Particularly, bulk viscosity obtained from the 
Gubser, Pufu and Rocha (GPR) formula for the GPR model [22] and the Improved Holo- 
graphic QCD model [23] coincides with the EO formula. The essential feature of the GPR 
formula (extracted from the holographic Kubo formula) is that it is suitable only for the 
models with one gravitational scalar field acting as the new radial coordinate. The exotic 



black hole model is the example of the phenomenological model with several gravitational 
scalar fields. 

We checked the validity of the Eling-Oz formula analytically for the exotic black holes 
in the high-temperature (conformal) limit. The formula is correct for the intermediate 
temperatures, the vicinity of the phase transition and for the temperatures up to ^ < 
2.75, where m is the mass associated with breaking of the conformal symmetry. ^ The 
correctness of the formula for exotic black holes extends the number of models for which 
the EO formula is valid for all energy scales. It would be interesting to explain this kind 
of universality of (1.2) (see also footnote [4] in the text). While preparing this manuscript, 

1 learned that the authors of [1] proved such a universality of the transport coefficients of 
the holographic plasmas [24]. They showed that the transport coefficients depend on the 
boundary conditions, but they are independent of the RG running from UV to IR. 

2 Bulk viscosity for the exotic hairy black holes 

The exotic black hole model is defined by the following effective (3+l)-dimensional gravi- 
tational action [2, 25]: 
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where g' is a coupling constant.^ Note that (j) induces a relevant deformation of the dual 
CFT by an operator Or and x is associated with an irrelevant operator Oj in the dual 
gauge theory. The last term in (2.1) involves mixing of Oi with Or under RG dynamics. 
The central charge of the UV fixed point is defined as [2] 

c=^, (2.2) 

This central charge should be understood as a measure of the degrees of freedom in the 
CFT, which is defined thermodynamically or via two-point correlation functions [26]. We 
demand the solution to be AdS^ asymptotically with translationary invariant horizon. For 
this purpose only the normalizable mode of Oi is nonzero near the boundary. 
The background geometry is defined as 

dsl = -ci{rf dt^+C2{r)^ [dxl+dxl]+C3irf dr\ ({> = 4>{r) , X = x(r) , (2.3) 

where r — )• oo corresponds to the AdS boundary. Then one can introduce a new radial 
coordinate x as follows 

l-x.SM, (2.4) 

C2[r) 

so that X —7- corresponds to the AdS boundary, and y = 1 — x — ?• corresponds to a 
horizon asymptotic. Afterwards, we introduce a{x) as 



^Original calculations of the bulk viscosity were done in [4, 5). 
^In numerical analysis we set g — —100. 



The equations of motion (EOMs) obtained from 2.1, with the background ansatz (2.3), 
define the following expansion of the model parameters 
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<^ =p^ X^^ + P2 x2/3 + —pfx + 0(x^/3) ^ (2.6) 

near the boundary x — t- 0+, and 

a = a(^a^ + a1y^ + 0(/)) , 'A = pS + O(y') , X = c|} + 0{y^) , (2.7) 

near the horizon. Up to the overall scaling factor a the thermodynamics of the black 
branes can be uniquely specified with 3 UV coefficients {pi,P2,Xi} ^^'^ 4 IR coefficients 

We use the integral of motion [25] 
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which arises after integration of the EOMs, to find the temperature T and the entropy 
density s of the black brane solution (2.3): 

- ^" (2.9) 
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s=^^ s = 4.W {4f , (2.10) 

In the dual picture pi can be interpreted as the coupling of the operator Or, P2 as the 
expectation value of Or and X4 as < Oj > (see [27] for the argumentation). Without the 
loss of generality, we choose the model with dim[C'r] = 2. Which means the combination 
pia should be fixed. 

For a given set of {a,pi} there is a discrete set of the remaining parameters 

{4, pIc'o} 

allowing us to characterize the thermodynamics of black branes suitable for Eling-Oz for- 
mula. One solution with Cq = describes the black brane without the condensate of the x 
field. All the other solutions have Cg 7^ and describe the "exotic black branes" [2]. This 
model is interesting because the transition occurs at the high temperatures (rather than the 
low temperatures). That is the irrelevant operator Oi obtains nonzero vacuum expectation 
value for T > Tc , spontaneously breaking a discrete Z2 symmetry of the model. Also, it 
was shown in [5] that all the exotic black branes contain a tachyonic quasinormal mode. 
Thus, they are dynamically unstable but thermodynamically stable, thereby, violating the 
correlated stability conjecture [28-30]. 

For the exotic black hole model formula (1.2) takes the following simple form 

^-l^ = ctr'\(^] +[^] I, (2.11) 
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Figure 1. (Colour online) Comparison of a scalar field at the horizon and the bulk viscosity 
computed using quasinormal modes with the high-temperature limit [2]. The dashed green lines 
represent the conformal limit. 



where Cs is a speed of sound defined by 



d(lnr) 
d{\xis)' 



(2.12) 



^^. Further, we will 



and r is an inversed dimensionless temperature, e.g. r = ^ or r 
check the validity of the formula for the different temperature regimes. 

2.1 Explicit analytical check of (1.2) in the conformal limit of a symmetric 
phase 

In this section we briefly repeat the main results for the thermodynamics of the model in 
the high temperature limit as presented in [25]. These results can be readily applied to 
supply evidence of (2.11) in the conformal limit. We demonstrate this below. For the sake 
of simplicity, we consider the symmetric phase only {x = 0). If we introduce the small 
deformation parameter 5 such that 
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where m is the mass associated with the deformation of CFT. Then one can expand the 
scalar field in the EOMs to leading order in 5 as ^ 



(t){y) = 6cj){y). 



(2.14) 



Afterwards, we solve the EOMs demanding regularity of the scalar field at the horizon and 
the first law of thermodynamics to fix the integration constants. ^ 
The expansion of the scalar field (j) near horizon (2.7) assumes that 
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^It is a natural small parameter to expand the solutions of EOMs defining the values of the scalar fields. 

*In general, p2 — the expectation value of Or is connected to the horizon data through the boundary 
conditions, regularity of the solution of the EOMs at the horizon and boundary. Particularly, in the 
conformal limit p2 cc Po- This can be the qualitative argument that the EO formula captures the UV data 
from the boundary. 




Figure 2. (Colour online) Comparison of the EO formula for exotic black holes with the computa- 
tions using quasinormal modes [4] far from criticality (the left plot) and in the vicinity of the critical 
point (the right plot). Red circles are for the "ordered" (symmetric) phase, purple diamonds are 
for the "disordered" (broken) phase. The dashed vertical green line represents the critical point of 
the theory T = T^. 



Eventually, it leads to the following expression for the speed of sound in the conformal 
limit 

cl = \-^5^ + 0{5'). (2.16) 

The results for the bulk viscosity in the conformal limit were discussed in [4]. For a 
symmetric phase at the high temperatures we have 
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A a result, the bulk viscosity should be proportional to the square of the deformation 
parameter, i.e., 

^ = -/■ (2.18) 

rj 4 

If we substitute (2.15) and (2.16) into (2.11) we recover exactly the same relation {t = 6 
in this case). This justifies the validity of (2.11) in the conformal limit. 

One can see that in the conformal limit the value of the scalar field at the horizon is 
proportional to ^. In Fig. 1 we use a linear approximation for the scalar field at the 
horizon to establish the connection between 5 and ^. In addition, on the same figure we 
compare the bulk viscosity evaluated from the sound waves attenuation coefficient with 
analytical result (2.18) in the conformal limit. 

2.2 Comparison of (1.2) aw^ay from criticality for symmetric and unstable 
phases 

Now it is possible to check the formula for an intermediate temperature regime. The 
original calculations of the thermodynamics and the bulk viscosity were done in [2, 4]. We 
use a cubic spline approximation for the values of the scalar fields at the horizon to get 
their derivatives with respect to the inversed temperature. Then we compare the validity 
of the Eling-Oz formula with respect to two arguments 5 and ^. 

Fig. 2 illustrates the absolute value of relative error between the bulk viscosity obtained 
from quasinormal mode method and the Eling-Oz formula. The error is slightly big for large 
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Figure 3. (Colour online) The square of a scalar field at the horizon Cg as function of the reduced 
temperature ^. The dashed green line is a linear fit to (cg )^ The right plot shows the bulk viscosity 
dependence on Cg. The dashed green line represents a linear fit to the log- log data. 



temperatures due to small values of the scalar fields (which increases numerical errors). But 
we have the analytic results from the previous section for this region of temperatures. The 

agreement is also worse in the vicinity of the critical point (due to large values of -37-)- In 
the next subsection we will improve the agreement in the critical regime. 

2.3 Comparison of (1.2) at criticality 

The critical behavior of the exotic black branes was discussed in [4]. The peculiar thing 
in the model is that the bulk-to-shear viscosity diverges in the broken phase at criticality. 
One can use more detailed data to check the formula close to criticality as it is done in the 
Fig. 2. Alternatively, we can proceed the semi- numerical analysis. 

In [4] the authors constructed an exotic model of the second order transition in d = 3 at 
finite temperature and zero chemical potentials. The corresponding conformal field theory 
in 2 + 1 dimensions is deformed by a relevant operator Or- The expectation value of Oi 
acts as the order parameter of the phase transition and it scales as |tp^ in the vicinity of 

the critical point, where t = ^ " . In Fig. 3 we plot In ( - I versus In (cq) with a dashed 
green line which fits the data 

2/ = -1.999(8)x- 9.427(8). (2.19) 

It is clear from the Fig. 3 that parameter cq has the same critical exponent as < Oj > 

Co oc |t|^/^ (2.20) 

For the speed of sound and bulk viscosity it was shown that 
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Whereas, formula (1.2) suggests that at criticality 



ctT^ 



EO 



Co 



T-T, 



oc \t\ 



(2.21) 



(2.22) 



which confirms the correctness of the EO formula at criticality. Also, we expect that the 
formula will be valid for other symmetry-broken phases. 
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Figure 4. Comparison of the EO formula for symmetric pliase witli tlie computations using 
quasinormal modes for the low temperatures. 

2.4 Validity of (1.2) for the low temperatures 

Let us proceed to the check of the Eling-Oz formula in the case when the stable phase 
is driven into the low temperature regime. First note that data for the thermodynamic 
parameters is more discrete, while the number of points for viscosity is significantly reduced 
(100 points altogether). Furthermore, Fig. 4 presents a plot of the relative error with 
respect to 6 and if. One can readily see that formula (2.11) is valid up to (5 < 2.75. 
Therefore, combining the results of [3], we expect that the formula is valid for the whole 
range of the temperatures. 
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